Pseudogap in the chain states of YBCO 
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As established by scanning tunneling microscopy (STM) cleaved surfaces of the high temperature super- 
conductor YBa2Cu207_5 develop charge density wave (CDW) modulations in the one-dimensional (ID) CuO 
chains. At the same time, no signatures of the CDW have been reported in the spectral function of the chain 
band previously studied by photoemission. We use soft X-ray angle resolved photoemission (SX-ARPES) 
to detect a chain-derived surface band that had not been detected in previous work. The 2/cp for the new 
surface band is found to be 0.55 A~\ which matches the wave vector of the CDW observed in direct space by 
STM. This reveals the relevance of the Fermi surface nesting for the formation of CDWs in the CuO chains in 
YBa2Cu207_5. In agreement with the short range nature of the CDW order the newly detected surface band 
exhibits a pseudogap, whose energy scale also corresponds to that observed by STM. 

PACS numbers: 71.45.Lr , 74.72.Kf, 79.60.-i, 79.60.Jv 



The CuO chain structure in YBa2Cu207_5 (Y-123) pro- 
vides a physical realization of a quasi-one-dimensional elec- 
tronic system w^ith non-vanishing coupling to the CUO2 
bilayers.^'^ One of the reasons w^hy ID electronic systems 
remain in the focus of solid state research is that even weak 
interactions transform the quasi-particles of Fermi liquid 
theory into collective excitations of density w^ave type.^~^ 
In particular, the peculiar topological structure of the chain 
Fermi surface (FS) makes this electronic system prone to 
formation of charge density wave via the Peierls instability.^ 
Indeed, charge density modulations along the Cu-0 chains 
have been extensively studied probing cleaved surfaces of Y- 
123 in direct space with scanning tunneling microscopy. 
It has been shown that Y-123 crystals cleave between the 
CuO chains and the BaO layer, so that the chains turn out 
to be the nearest to the surface building block.^'^'^^'^^ Both 
the earliest^ and the most recent studies^"^ present a con- 
sistent picture of a CDW, appearing as corrugations of the 
electronic density with a short correlation range of about 
40 A and a period between 9 and 14 A, depending on the 
sample stoichiometry. 

In the case of 2D systems, the effect of CDW on the elec- 
tronic states in the reciprocal space has been examined in 
great detail with modern angle resolved spectroscopy, which 
provides both momentum and energy resolution when mea- 
suring the one particle spectral function. Occurrences of 
the pseudogap for an incommensurate (or short range or- 
dered) state, which finally develops into a true CDW band 
gap below T^dw^ ^re well documented. In the case of 
ID systems other than Y-123, modifications to the spectral 
function with the onset of the CDW state have been detected 
as well.^^~^^ the emergence of a pseudogap-like state is also 
expected in theory and can be understood as a result of a 
fluctuating Peierls order parameter.^^"^^ 

At the same time, despite the abundant evidence for 
charge modulations in direct space provided by STM, 
Y-123 chains seem to exhibit neither a pseudogap nor 
the Tomonaga-Luttinger behavior observed in the related 



PrBa2Cu408.^^'^^ Momentum resolved spectra measured 
from as-cleaved Y-123 surfaces,^^"^^ as well as those mea- 
sured on the in situ doped ones"^^ manifest neither folding 
nor notable suppression of the spectral weight at the Fermi 
level. That is, the spectral features that could be regarded as 
signatures of the CDW state are absent in spectra of Y-123. 
Therefore, it remains unclear why the two complementary 
methods (STM and ARPES) deliver such conflicting results. 

To address this problem we have investigated high qual- 
ity Y-123 crystals using modern SX-ARPES. the photoemis- 
sion data for this study were collected at the recently built 
high-resolution soft X-ray beamline ADRESS at the Swiss 
Light Source. This beamline delivers photon flux up to 
10^^ ph/s in an energy window of 0.01% of photon energy. 
Such a high flux allows one to break through the problem 
of dramatic reduction of the valence band cross section 
at high photon energies, the samples were mounted on a 
low-temperature goniometric manipulator (CARVING) with 
three angular degrees of freedom and cleaved in situ in ul- 
trahigh vacuum with base pressure better than 5-10"^^ mBar 
at r=10K, the same temperature at which all the spectra 
were acquired. The energy resolution depends on the kinetic 
energy and will be stated in the figure captions, details of 
the SX-ARPES station will be published elsewhere.^^ The 
experiments were performed using high-quality single crys- 
tals of Ba2Cu305 5 (T^^61 K) from the same batch as in the 
recent INS study. The crystals were synthesized by the 
solution-growth technique, annealed to the desired oxygen 
doping, and detwinned by applying uniaxial mechanical 
stress at elevated temperature.^^ 

As a result of the 'surface polar catastrophe'^"^ the elec- 
tronic structure of the Y-123 surface layer is known to be 
different from that of the bulk.^^ Generally, there are two 
major approaches that have been used in photoemission 
practice to highlight bulk states against the surface ones. 
Surface aging is the first option, which basically relies on 
the destruction of surface sates. ^^'^^ The other approach 
employs the variation of photon excitation energy and po- 
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FIG. 1: Sample aging, (a), (b) FS maps measured directly after 
cleavage and 28 hours later (c), (d) Corresponding fresh and aged 
intensity distributions for constant energy cuts at £ = £p — 100 meV 
(e), (f) Energy-momentum intensity distributions along the F-S 
direction. The white squares are the result of MDC fits. Dark 
corresponds to high photoemission intensity. Spectra measured at 
r = 10 K with 38 meV energy resolution using linearly s -polarised 
light (i.e. no out of plane component of the polarization vector). 



larization, v\Ahich affect the inelastic escape depth of exited 
electrons"^^""^^ and the ratio of photoemission matrix ele- 
ments for bulk and surface states. '^^''^'^ Also the rapid varia- 
tion in the elastic scattering rate v\^ith the kinetic energy of 
the excited photoelectrons may play a decisive role."^"^'"^^ 

We start out the discussion of experimental data with an 
illustration of the first option. Fig. 1 contains a comparison 
of the Fermi surfaces measured from the freshly cleaved 
sample to the same measurement done 28 hours after the 
cleavage. The freshly cleaved surface results in a typical pic- 
ture of bilayer split FS contours corresponding to the CUO2 
plane states (rounded double squares centered at S point) 
and the quasi-one-dimensional chain band (faint features 
running parallel to the axis) . After surface degradation 
there are some notable changes. 

First, the overdoped CUO2 plane bands are replaced v\^ith 
the 'Fermi arc' features similar to those observed in refs. 46, 
47. The appearance of the 'Fermi arcs' is also accompanied 
by a /Cp shift (Fig. l(e-f)). For the bonding band, it can 



be estimated from the MDC fits, and amounts to about 
5k^ ^ 0.04 A along the F-S direction. Assuming the same 
shrinking over the whole Fermi surface and an average FS 
radius ^ 0.55 A, one may estimate the decrease in hole 
doping as 5p ^ 4nk^5k^/S^^ ^ 0.1, which brings us from the 
typically overdoped regime with p ^ 0.3 (refs. 30,31,38,48) 
under the superconducting dome with p ^ 0.2. Note also 
the so called 'waterfalls' and a high energy kink at about 
250 mey which have recently generated an avalanche of 
publications. "^^"^^ These issues are beyond the scope of the 
current study, and we do not further elaborate on them 
here. 

Second, the relative intensity of the ID chain band as 
compared to the 2D bands is substantially increased, while 
the distance between the chain Fermi crossings, 2A:p, re- 
mains practically the same. This suggests that the observed 
chain band is most likely a 'bulk' feature, in a sense that it 
originates from the chain structure protected by at least one 
CUO2 bilayer, and not from the cleaved chains. Since the 
upper neighboring CUO2 bilayers turn out to be overdoped, 
the electronic structure of these chains still differs from LDA 
band structure calculation,^ so in the following we will refer 
to these chains as subsurface ones, in order to contrast this 
structure with the chains through which the cleavage takes 
place and the true bulk chains. 

This observation is consistent with the interpretation 
given in ref. 32, which was based on the pattern of cir- 
cular dichroism. This work has shown that the chain sig- 
nal observed in ARPES experiment mainly originates from 
the nearest to the surface undisturbed chains protected by 
the overdoped CUO2 bilayer, so that the surface aging is 
expected to enhance the photoemission from these chain 
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FIG. 2: Excitation energy dependence for the cut passing through 
the r point parallel to the ky axis. The left panel shows the inten- 
sity integrated around the Fermi level, E = Ep. The right image 
contains intensity at 200 meV below the Fermi level, demonstrat- 
ing another set of chain features for hv <920 eV Spectra were 
measured using p -polarized Hght with overall energy resolution of 
125 meV 
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states as compared to other surface features. Though there 
are intensity variations along the chain band, it is notewor- 
thy that the band exhibits neither pseudogap nor folding 
phenomena. As can be seen in Fig. l(a-d), there are no 
indications for a CDW state for the freshly cleaved and aged 
surfaces. 

Now we turn to the discussion of the second option. 
Searching for the signatures of a CDW in Y-123 spectra 
we have measured the excitation energy dependence for 
the Y-F-Y cut, the one where only a parabolic chain band 
is expected to cross the FL.^ Fig. 2 contains the intensity 
distribution for that direction plotted at the FL (a) and 
200 meV below the FL (b) . While the data in panel (a) ex- 
hibit only two parallel features, corresponding to the di/cp 
crossings of the two branches of the parabolic chain band, 
the intensity distribution in Fig. 2(b) reveals another set of 
features running parallel with a slightly larger separation. 
The new features are well visible at low excitation energies 
hv <920ey therefore based on the universal dependence 
of the escape depth one may assume the new feature to be 
a surface related one. However, the decrease of the photoe- 
mission intensity from the surface chain band as compared 
to the subsurface ones occurs rather abruptly, whereas the 

1 /2 

universal escape depth curve predicts smooth oc E^!^ de- 
pendence, where E^^^ is the kinetic energy of the excited 
photoelectrons. Therefore this observation cannot be solely 
attributed to the electron escape depth, so also the hv de- 
pendence of the photoemission matrix elements needs to be 
taken into account. In contrasting surface photoemission 
with bulk emission, particular attention has to be payed 
to the surface induced light fields, which are likely to be 
responsible for the different hv dependence of the matrix 
elements for surface and bulk localized states. For in- 
stance, the absence of the new surface chain states in the 
spectra shown in Fig. 1 would be consistent with the fact 
that these spectra were measured with s -polarized light. For 
this polarization there is no component of the exciting field 
A(r) perpendicular to the sample surface, and hence the 
surface related term in the photoemission matrix element 
(f|AV + I divA|i) is negligible, rendering the surface emis- 
sion effects irrelevant as compared to the case of p -polarized 
light. 

To establish the dimensionality and topology of this pre- 
viously overlooked surface state, in Fig. 3 (a,b) we plot a 
FS measured with 620 eV photons and corresponding iso- 
energy intensity distribution taken 100 meV below the FL. 
According to the hv dependence, the new feature is expected 
to have significant contribution to the experimental spectra 
at this energy. As can be seen, the new band practically does 
not contribute spectral weight at the FL (pseudogapped), 
though at 100 meV binding energy the band gives rise to two 
ID traces similar to those spawned by the chain band. This 
suggests that the subsurface chains and the pseudogapped 
surface feature stem from homologous bands with, probably, 
different spatial localization. Indeed, a detailed analysis of 
the LDA band structure^ shows that there are no other ID 
bands except for the chains that the pseudogapped feature 
could be attributed to. Further, dispersion cannot be 
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FIG. 3: (a), (b) Momentum intensity distribution at the Fermi 
level and at 0.1 eV below the FL. The right image reveals another 
ID band, which contributes practically no spectral weight at the 
Fermi level. Spectra measured with p-polarized light and energy 
resolution of 90 meV 



responsible for the appearance of the second pseudogapped 
chain band in the spectrum, since along the Y-F-Y direction 
the calculated variation with (±4% k^) is notably 
smaller than the observed splitting. 

In Fig. 4 (a, b) we plot the energy-momentum intensity 
distribution for the Y-F-Y direction. Panel (b) contains the 
same data set, but each MDC making up the image has 
been normalized to a fixed value, in order to show the band 
dispersion in the region that is 'overexposed' in panel (a) . 

The inner band with 2k^ ~0.28 A~^ is the chain band that 
we attribute here to the subsurface chain states, and the 
outer one, with 2k^ ~0.55 A,~^ is the newly observed band, 
which we believe to be a surface localized chain band. The 
suppression of the spectral weight down to about 250 meV 
for the surface chains as compared to the bulk ones is clearly 
visible. In the panel (c) we additionally compare two EDC's 
taken at the effective k^ of the bulk and surface bands, which 
once again demonstrate the pseudogap-like suppression of 
the spectral weight for the surface band. While for the sub- 
surface chains the spectral intensity grows approximately 
as the Fermi function convoluted with experimental resolu- 
tion, for the surface band a gradual growth of intensity is 
observed. The estimated energy scale is about 200-300 meV 

In the case of true CDW gap, using model calculations, it 
has been demonstrated that new states arise on the order 
of Ep di A. For these energies, the CDW component of the 
local DOS was shown to have a characteristic amplitude 
and phase. Extending this results on the pseudogap, we 
note that the pseudogap nicely agrees with the STM results 
of ref. 7 and ref. 12. Namely, the charge corrugations are 
visible for bias voltages up to 310 meV and practically vanish 
for bias voltages higher than 480 meV In this regard we can 
also conclude that another small gap (~ 20 meV) observed 
in STM spectra and attributed either to CDW gap or to 
proximity induced superconducting gap,^ is likely to be a 
superconductivity related one. 
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FIG. 4: (a) Energy-momentum cut along 
the Y-r-Y direction demonstrating two 
types of chain bands, (b) The same as 
previous, but each MDC in the image 
has been normaHzed to one. The small 
square symbols denote the peak positions 
obtained in a four-Lorentzian MDC fit. 
The insent shows typical fit for the MDC 
at £B=200meV The obtained experimen- 
tal dispersions are fit to parabolic bands 
(solid white lines), (c) /cp EDC's for the 
inner chain band (black line) and outer 
chain band (grey line) . Polarization and 
resolution are the same as in Fig. 3. 



It is also informative to see how the observed band disper- 
sion compares to the COW modulations measured in STM ex- 
periment. Simple theory predicts that the COW wave vector 
should be twice the Fermi vector, q^^y^ = 27i/Acdw = 2A:p. 
To estimate from the photoemission spectrum where the 
pseudogapped chain band would cross the FL, if there were 
no CDW instability, we use the result of an MDC fit per- 
formed for energies below the gapped region (white symbols 
in Fig. 4(b) and extrapolate the experimental points up to 
the Fermi level assuming a parabolic dispersion. The value 
thus obtained is 2k^ = 0.55 A~^, which yields ?icDw = 11.4 A. 
This value nicely compares to the X= 11.2 A, observed in 
Zn substituted Y-123.^^ It is also within the range of CDW 
periods from 9 to 14 A reported for a series of YBa2Cu207_5 
samples with varying stoichiometry/^ approximately corre- 
sponding to the sample with 5 = 0.35 used in that study. 

At the same time for the subsurface chains with 2k^ = 
0.28 A"\ so the expected CDW period would be about 22 A, 
which obviously does not mach the STM data. Together 
with the absence of the (pseudo)gap this once again sug- 
gests that CDW modulations, which have been studied in 
so many details by STM, are due to the 2k^ = 0.55 A~^ 
instability occurring in the previously overlooked surface 
chain band. Disorder in the chain structure results in a short 
range character of the CDW and the development of the 
pseudogap in the chain band. 

It is remarkable that the chains at the surface develop 
the CDW state, while the neighboring chains sandwiched 
between the overdoped and bulk CUO2 layers display no 
notable signatures of the CDW. One reason for this could 
be a trivial difference in phonon modes available at the 
surface and in the bulk. However commensuration effects 



appear to be of no less importance. As compared to the 
incommensurate state, there is an additional energy low- 
ering associated with the commensurability of the CDW 
state, which is not taken into account in the simplest vari- 
ant of the CDW theory. ^'^^ The energy gain is given by 

= Here D is the band width, A 

is the CDW gap, is the electron-phonon coupling 
constant, n(e^) - density of states at the Fermi level, and 
M = AcDw/ a is the commensurability factor As can be seen 
the correction is most significant at low M values. Con- 
sidering this correction, for the surface chains we have 
M ^ 11.4/3.89 ^ 2.94, which is indistinguishable from 
a commensurate modulation with M = 3. For the subsur- 
face chains M = 5.77 is quite large and differs significantly 
from the nearest commensurate value M = 6. Consequently 
the E^^^^ correction is negligible in this case, which obvi- 
ously should affect the energy balance, making the CDW 
phase less favorable for the subsurface chains. 

In conclusion, we have identified the origin of the CDW 
observed at the surface of Y-123 by STM by detecting a pre- 
viously overlooked ID surface band in Y-123, which brings 
into agreement two complementary methods: ARPES and 
STM. As a result of a Peierls instability, the surface CuO 
chain band develops a short range CDW state, which results 
in the appearance of a pseudogap in the one particle exci- 
tation spectrum detected in our ARPES experiment. Both 
energy and momentum scales measured in STM and ARPES 
are found to be in agreement with each other 

We thank F. Dubi and C. Hess for technical support during 
the measurements. The project was supported in part by 
DFG grant ZA 654/1-1. 
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